Inhomogeneous Diastereomeric Composition of Mongersen Antisense Phosphorothioate Oligonucleotide Preparations and Related Pharmacological Activity Impairment

Mongersen is a 21-mer antisense oligonucleotide designed to downregulate Mothers against decapentaplegic homolog 7 (SMAD7) expression to treat Crohn's disease. Mongersen was manufactured in numerous batches at different scales during several years of clinical development, which all appeared identical, using common physicochemical analytical techniques, while only phosphorous-31 nuclear magnetic resonance (P-31-NMR) in solution showed marked differences. Close-up analysis of 27 mongersen batches revealed marked differences in SMAD7 downregulation in a cell-based assay. Principal component analysis of P-31-NMR profiles showed strong correlation with SMAD7 downregulation and, therefore, with pharmacological efficacy in vitro. Mongersen contains 20 phosphorothioate (PS) linkages, whose chirality (Rp/Sp) was not controlled during manufacturing. A different diastereomeric composition throughout batches would lead to superimposable analytical data, but to distinct P-31-NMR profiles, as indeed we found. We tentatively suggest that this may be the origin of different biological activity. As similar manifolds are expected for other PS-based oligonucleotides, the protocol described here provides a general method to identify PS chirality issues and a chemometric tool to score each preparation for this elusive feature

Metropolis simulations of Met-Enkephalin with solvent-accessible area parameterizations

We investigate the solvent-accessible area method by means of Metropolis simulations of the brain peptide Met-Enkephalin at 300$K$. For the energy function ECEPP/2 nine atomic solvation parameter (ASP) sets are studied. The simulations are compared with one another, with simulations with a distance dependent electrostatic permittivity $\epsilon (r)$, and with vacuum simulations ($\epsilon =2$). Parallel tempering and the biased Metropolis techniques RM$_1$ are employed and their performance is evaluated. The measured observables include energy and dihedral probability densities (pds), integrated autocorrelation times, and acceptance rates. Two of the ASP sets turn out to be unsuitable for these simulations. For all other systems selected configurations are minimized in search of the global energy minima, which are found for the vacuum and the $\epsilon(r)$ system, but for none of the ASP models. Other observables show a remarkable dependence on the ASPs. In particular, we find three ASP sets for which the autocorrelations at 300$$K are considerably smaller than for vacuum simulations.Comment: 10 pages and 8 figure

Water activity in liquid food systems : A molecular scale interpretation

Water activity has historically been and continues to be recognised as a key concept in the area of food science. Despite its ubiquitous utilisation, it still appears as though there is confusion concerning its molecular basis, even within simple, single component solutions. Here, by close examination of the well-known Norrish equation and subsequent application of a rigorous statistical theory, we are able to shed light on such an origin. Our findings highlight the importance of solute-solute interactions thus questioning traditional, empirically based “free water” and “water structure” hypotheses. Conversely, they support the theory of “solute hydration and clustering” which advocates the interplay of solute-solute and solute-water interactions but crucially, they do so in a manner which is free of any estimations and approximations

Effect of the integration method on the accuracy and computational efficiency of free energy calculations using thermodynamic integration

Although calculations of free energy using molecular dynamics simulations have gained significant importance in the chemical and biochemical fields, they still remain quite computationally intensive. Furthermore, when using thermodynamic integration, numerical evaluation of the integral of the Hamiltonian with respect to the coupling parameter may introduce unwanted errors in the free energy. In this paper, we compare the performance of two numerical integration techniques-the trapezoidal and Simpson's rules and propose a new method, based on the analytic integration of physically based fitting functions that are able to accurately describe the behavior of the data. We develop and test our methodology by performing detailed studies on two prototype systems, hydrated methane and hydrated methanol, and treat Lennard-Jones and electrostatic contributions separately. We conclude that the widely used trapezoidal rule may introduce systematic errors in the calculation, but these errors are reduced if Simpson's rule is employed, at least for the electrostatic component. Furthermore, by fitting thermodynamic integration data, we are able to obtain precise free energy estimates using significantly fewer data points (5 intermediate states for the electrostatic component and 11 for the Lennard-Jones term), thus significantly decreasing the associated computational cost. Our method and improved protocol were successfully validated by computing the free energy of more complex systems hydration of 2-methylbutanol and of 4-nitrophenol-thus paving the way for widespread use in solvation free energy calculations of drug molecules